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ABSTRACT 

Using a generalized non-spherical, multi-transonic accretion flow model, we analytically calculate the normalized 
QPO frequency V qpo of galactic black hole candidates in terms of dynamical flow variables and self-consistently 
study the dependence of v qpo on such variables. Our results are in fairly close agreement with the observed 
QPO frequencies of GRS 1915+105. We find that V qpo is quite sensitive to various parameters describing the 
black hole accretion flow containing dissipative and non-dissipative shock waves. Thus the QPO phenomena is, 
indeed, regulated by 'shocked' black hole accretion, and, for the first time, we establish a definitive connection 
between the QPO frequency and the properties of advective BH accretion flows. This information may provide 
the explanation of some important observations of galactic micro quasars. 

Subject headings: accretion, accretion disks — black hole physics — hydrodynamics — shock waves — QPO 



1. INTRODUCTION 

It has been established in recent years that a number of galac- 
tic black hole (BH) candidates show quasi-periodic oscillation 
(QPO) behaviours spanning a wide range of frequency starting 
from a few tens of milliHertz to a few hundred Hz (see Cui 1999 
for a review). The characteristic features of these QPOs sug- 
gest that such oscillations are a diagnostic of the accretion pro- 
cesses in the inner regions of the accretion discs around the cen- 
tral compact objects, like shock oscillations in sub-Keplerian 
advective accretion discs (Chakrabarti & Titarchuk 1995, CT 
hereafter, Molteni, Sponholz & Chakrabarti 1996, MSC here- 
after, Kazanas, Hua & Titarchuk 1997, Titarchuk, Lapidus & 
Muslimov 1998, Chakrabarti & Manickam 2000, CM here- 
after). The galactic microquasar GRS 1915+105 shows several 
properties which indicate that QPOs may be generated because 
of the oscillation of the hot, dense inner regions of black hole 
accretion discs (Paul et al. 1998, Muno, Morgan & Remillard 
1999, Yadav et al. 1999, Rutledge et al. 1999, Smith, Heindl & 
Swank 2002). Numerical simulation of MSC revealed the fact 
that for shocked non-spherical advective BH accretion flow, if 
the inflow and cooling parameters are properly tuned, the post- 
shock cooling time becomes comparable to the advection time- 
scale from the shock. A shocked flow may continuously try to 
adjust itself, resulting in the quasi-periodic radial oscillation of 
the shock around its mean steady-state solution, with an oscilla- 
tion frequency roughly equal to the reciprocal of the post-shock 
advection time scale r a </,.. Such oscillations do not depend on 
the details of the cooling mechanism (CM) and may cause the 
quasi-periodic variation of the luminosity emerging out of the 
post-shock region of the accretion disc by dynamically modu- 
lating the size and the radiative properties of the Comptonizing 
cloud (CT). One thus expects that the shock formation in multi- 
transonic advective BH accretion discs may play a crucial role 
in producing QPOs. 

In this letter, we study how sensitive shock oscillation is to 
the parameters governing the BH accretion, which has never 



been done in any of the works existing in the literature. Our 
main motivation in this letter is to provide a generalized self- 
consistent computation of various shock parameters V s h in 
terms of fundamental inflow variables V ace for accretion in 
all available post-Newtonian pseudo-Schwarzschild BH poten- 
tials, and then to use these parameters to estimate the value of 
Vqpo- We also study the dependence of v qpo on all V acc and V s h- 
We assume that the shock solutions we obtain in this letter, sat- 
isfy the criteria of fine tuning inbetween the post-shock advec- 
tion and cooling time scales so that shocked flow oscillates with 
a frequency ~ 1 /r a d v even for realistic flows with all dissipative 
terms included. Our generalized formalism assures that our re- 
sults are not just an artifact of a particular type of potential only 
and inclusion of every BH potential allows a substantially ex- 
tended zone of parameter space allowing for the possibility of 
shock induced QPOs. Also, there are possibilities that in future 
someone may come up with a pseudo-Schwarzschild potential 
better than any of the BH potentials available in the literature, 
which will be the best approximation for complete general rel- 
ativistic investigation of multi-transonic shocked flow. In such 
case, if one already formulates a generalized model for multi- 
transonic shocked accretion disc for any arbitrary potential, ex- 
actly what we have done in this paper, then that generalized 
model will be able to readily accommodate that new poten- 
tial without requiring any significant change in the fundamental 
structure of the formulation and solution scheme of the model 
and we need not worry about providing any new scheme exclu- 
sively valid only for that new potential, if any. 

2. GOVERNING EQUATIONS, SOLUTION PROCEDURE AND 
RESULTS 

Following are the four pseudo-potentials on which we con- 
centrate in this work: 
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In this work, any quantity changing discontinuously at shock location is termed as shock parameters. 
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where r is the radial co-ordinate scaled in units of 
Schwarzschild radius r g (=2GM BH / 'c 2 ). Throughout this paper 
the radial distances and velocities are scaled in units of r g and c 
respectively and all other derived quantities are scaled accord- 
ingly. Also G = c = M BH = 1 is used. In eq. (1), $1 and $ 2 
had been proposed by Paczyhski & Witta (1980) and Nowak & 
Wagoner (1991) respectively, and $3 and $4 are suggested by 
Artemova et al (1996). For these potentials, Das (2002, D02 
hereafter) and Das, Pendharkar & Mitra (2003, DPM hereafter) 
have provided a unified approach for shock solutions in multi- 
transonic black hole accretion disc for non-dissipative shocks in 
polytropic flow (D02) as well as in dissipative shocks in isother- 
mal flows (DPM). We follow the formalism presented in D02 
and DPM to calculate V s h for evaluating v qpo . 

2. 1 . Polytropic accretion ; Rankine-Hugoniot shock 

For this case, V ace = {£, A, 7}, where £ is the specific energy 
of the flow, A is the specific flow angular momentum, and 7 is 
the adiabatic index of the flow. We calculate the shock loca- 
tion r s h, along with other shock parameters like shock strength 
Si(=M-/M+, M being the Mach number of the flow), shock 
compression ratio R comp (=p + / p-, p being the local matter den- 
sity of accretion disc, related to its vertically integrated value 
£ as p = T,/I n h(r), h(r) being the local disc height and /„ = 
(2"n!) 2 /(2«+l)!, n = 1/(7-1), see Matsumoto et al. 1984) 
and entropy enhancement rate at shock (3 (=A4+ /A4-, M. be- 
ing the entropy accretion rate) as a function of V ace - For low A 
flow with moderate or strong shocks, the post-shock velocity is 

proportional to /? com/ ,r^ 2 (see also CM) and one obtains: 

1 _ A 

Tad, (~ r sh /u sh ) R comp r2 Z 

where A is a scaling constant. Let II be any physical quantity 
which changes discontinuously at the shock. We define the ratio 
of the post to pre-shock value of IT to be 5Rn = jp- We multiply 
the numerator and denominator of eq. (2) by R CO m P and use the 
following relations: 
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to obtain the QPO frequency as: 
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for GRS 1915+105. Hereafter, we define the normalized QPO 
frequency as v qpo = v qpo j(fl. 

In figure 1., we plot V qpo as a function of various fundamental 
'Pace- While the v qpo is plotted along Z axis (denoted by v in 
the figure), specific energy £ and specific angular momentum A 
are plotted along X and Y axis respectively for 3>i(A), $2(6), 
$3(C) and $4(D). Each surface represents the plot for a fixed 7. 
The left most surface is drawn for 7 = 5 and successive surfaces 
to the right are drawn for (7 + n A7) where A7 = 0.1 and n is an 
integer. For^ as well as for $4, we obtain the figure for n= 1,2; 
while for $3 we use n=l and for $2 we use n=0. It is observed 
from the figure that V qpo anti-correlates with A and 7. This in- 
dicates that the ultra-relativistic flow 2 produces stronger QPO 
compared to the purely-non relativistic flows when other ini- 
tial parameters are unchanged. Weakly rotating flows are better 
candidates for producing high frequency oscillation; this is due 
to the fact that the shock location anti-correlates with V qpo but 
correlates with A (see Fig. 2). The lower the be the value of 
angular momentum, the closer will be the post shock region to 
the event horizon and more amount of the gravitational energy 
will be available to dump on to the oscillation energy in the 
post shock region. For most of the cases, v qpo correlates with 
£ while v qpo anti-correlates with £ for relatively fewer cases. 
Such aV qpo -£ (anti) correlation are supported by observational 
evidences that most QPOs strengthen towards higher energies 
(van der Klis 1995; Cui et al. 1997; Morgan et al. 1997; Remil- 
lard et al. 1998), although some follow the opposite trend (Mor- 
gan et al. 1997; see also Remillard et al. 1999 and Rutledge et 



al. 1999 for other possible cases). If ' 



where W and 9 are the vertically integrated fluid pressure and 
the flow temperature respectively. £f h is the ratio of flow kinetic 
energy to the total thermal energy content of the flow. All rel- 
evant II and 5ftn could be interpreted as V s h- For a particular 
set of V a cc for which shock forms for any specific <&„ all 5Rns 
could be computed in terms of V ace for all available <J>,s, hence 
we claim that the expression for v qpo in eq. (4) is more general 
compared to the QPO frequency obtained by any other previous 
work. A suitable choice of A provides the maximum value of 
v qpo (for any V acc in all $,) to be 67 Hz for one Eddington ac- 
cretion rate onto a 10M Q BH, which is the representative case 

2 By the term 'ultra-relativistic' and 'purely non-relativistic' we mean a flow with 7 
1992. 
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maximum possible value of normalized QPO frequency for a 
particular astrophysical source for any specific &i(r), we find 
that: 



(5) 



For ultra-relativistic shocked flow in certain range of £ and A, 
in Fig. 2, we plot v qpo (plotted along the Y axis and denoted by 
fnorm in the figure) as a function of various Vsh for all four $, 
marked in the figure. To plot in the same figure, the post shock 
temperature T + (dotted line) and the post shock vertically inte- 
grated density E + (long dashed line) is scaled as T+ — > a,T+ and 
£+ — ► where a, and b[ are scaling constants different for 
four different ®i(r). Similar figures could also be drawn for any 
range of {£ ,A} for which a shocked flow physically exists. It 
is evident from the figure that v qpo non-linearly anti-correlates 
with the shock location whereas it non linearly correlates with 
the post shock temperature and post shock density. This result 
may have important consequences in connection to the jet for- 
mation mechanism for galactic relativistic sources. It has been 
shown that (Das 1998, Das & Chakrabarti 1999) the hot and 
dense post shock region of black hole accretion disc may gener- 
ate accretion-powered galactic and extra-galactic outflows and 
the barionic load of such outflows increases with the increment 
of post shock flow temperature and density. From Fig. 2, we 
conclude that there must be some relationship inbetween the 
normalized QPO frequencies and the amount of baryonic con- 
tent of outflows from galactic microquasars; details of such in- 
vestigation is reported elsewhere (Das, Rao & Vadewale 2003). 
In very high states of the QPOs in LMXB, radiation drag might 

= i and 7=4 respectively, according to the terminology used in Frank et. al. 
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play a crucial role in accretion close to the event horizon and 
the QPO might be caused by the radiatively driven oscillations 
in a quasi-spherical accretion flow (Fortner et al. 1989). The 
Vq po — T + co-relation presented in Fig. 2 is thus expected to sup- 
port such QPO states because the quasi-spherical post-shock 
region in our model possesses more thermal energy as it forms 
closer to the event horizon. 

2.2. Isothermal accretion; dissipative shocks 

For such flow, we calculate all V s i t s and 3?^ in terms of 
V acc ={^,T}, where T is the constant flow temperature. The 
amount of energy dissipation at shock may become as much as 
five to fifteen (at most) percent of the rest mass energy of the 
accreting material and may be computed as: 



1 kT 
Ae=-(M 2 -Ml)-—lnft 1 

2 v " +; iim H 



(6) 



This huge amount of dissipated energy may be used to power 
the QPO and/or to produce and power the strong X-ray flares 
from the BH environment. We can calculate the normalized 
QPO frequency for this case as: 



exp 



Aefim H 
kT 
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In figure 3. we plot the variation of v qpo (plotted along the 
Z axis and is marked as v in the figure) with specific angu- 
lar momentum A (plotted along X axis) and flow temperature 
(plotted along Y in the unit of 7io (=T/10 10 Degree K) and is 
marked as r in the figure) for four pseudo potentials <J>i(A), 
$2(B), $3(0 and ( E>4(£>). It is observed from the figure that 
V qpo anti-correlates with A and correlates with r, and: 
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which is exactly same as that of polytropic flow, see eq. (5). 
This clearly indicates that shocked accretion flow in $3 gener- 
ates the most intensely oscillating QPOs. This indicates that the 
high frequency QPOs observed for any astrophysical sources 
may, perhaps, be better explained using $3, while ($1 , $4) and 
$2 are appropriate to explain the intermediate and low fre- 
quency QPOs, respectively. For shocked flow in certain range 
of A and T, figure 4 represents the variation of v qpo (plotted 
along the Y axis and is denoted by f noml in the figure) with 



shock location r s >, (solid line), shock compression ratio 5Rs 
(dotted line) and energy dissipation at shock Ae (long dashed 
line). Similar figures could also be drawn for any range of 
{A,T} for which a shocked flow physically exists. In the fig- 
ure, and Ae are scaled as 3?s — ► c,3?s and Ae — > <i,Ae 
where c, and <f, are scaling constants different for four different 
$i(r). It is observed that v qpo non-linearly anti-correlates with 
the shock location which indicates the generation of stronger 
oscillations closer to the event horizon as is observed in real- 
ity (Cui 1999). Also V qpo correlates with 5Rs and Ae which 
indicates that the stronger is the shock (as 5ft s = <S,-, 5, being 
the shock strength, see DPM), more intense is the oscillation 
and higher is the amount of dissipated energy at the shock loca- 
tions. Co-relation in between ~P qpo and Ae for isothermal flow 
may help in understanding the characteristic features of the du- 
ration of the quiescent state of the transient X ray source GRS 
1915+105 (Das & Rao, in preparation). 

3. CONCLUSION 

In this work we have examined the dependence of normal- 
ized QPO frequencies V qpo on various parameters describing 
the generalized multi-transonic BH accretion flow with dissi- 
pative and non-dissipative shock waves. We found that V qpo is 
quite sensitive on most of the parameters responsible for pro- 
ducing the post shock flow, the anticipated oscillation of which 
may resulting the QPOs. This indicates that the QPOs in galac- 
tic sources are inherently controlled by shock dominated ac- 
cretion flows. The fact that v qpo is sensitive to various flow 
parameters found in this work reinforces the belief that our cal- 
culation may be useful in future to understand various obser- 
vational features of QPO in galactic BH candidates, a few of 
which have already been addressed in §2.1 and §2.2. We per- 
formed our calculations in a generalized way by incorporating 
all available pseudo-Schwarzschild black hole potentials. How- 
ever, one limitation of the work presented here is that we could 
not postulate any observational method that might convincingly 
discriminate among our results obtained using the various po- 
tentials. 
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Figure Captions 



Fig. 1 : For polytropic accretion with RH shocks, varitaion of normalized QPO frequencies (indicated by v in the figure) with specific 
flow energy £ and angular momentum A for <&i(A), $2(6), ^(C) and $4(D). The figure is drawn for different 7 with the interval 
A 7 =0.1. 

Fig. 2: For polytropic accretion with RH shocks, varitaion of normalized QPO frequencies (indicated by f norm in the figure) with 
shock location r s h, post shock flow temperature T + and post shock integrated density £+. 

Fig. 3: For isothermal accretion with dissipative shocks, varitaion of normalized QPO frequencies (indicated by v in the figure) with 
specific flow temperature in units of 7io = T/IO 10 Degree K (indicated by r in the figure) and angular momentum A for <&i(A), $2(B), 
$ 3 (C) and $ 4 (D). 

Fig. 4: For isothermal accretion with dissipative shocks, varitaion of normalized QPO frequencies (indicated by f norm in the figure) 
with the shock location r s h, shock compression ratio ( = ^ ) an d the amount of energy dissipation at shock Ae. 
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